clear all
use "./meatall.dta"

egen trall = rowtotal(tr201 tr202 tr203 tr207 tr20610 tr20629 tr20630 tr20649)
gen lntrall = ln(trall)

* udb: CIF value in ud (US Dollar) 
gen udp = exp(lnud207) 

* tariff levied
gen udpt = udp*(1 + tr207/100)

* Poultry meat
gen trr207 = (1 + tr207/100)
gen trr203 = (1 + tr203/100)
gen trr201 = (1 + tr201/100)
gen trr202 = (1 + tr202/100)
gen lntrr207 = ln(trr207)
gen lntrr203 = ln(trr203)
gen lnxrjpy = ln(xrjpy)

* Instruments
gen lntxrp = lntrr207 + lnxrjpy //iv
gen lnudp = ln(udp) // Log CIF value in USD 
gen lnjypt = lnudp + lntxrp // Log expenditure incl tariff in JPY
gen jypt = exp(lnjypt)

* kgb: CIF in kg (quantity) 
gen kgp = exp(lnkg207)

* meat
gen lnkgp = ln(kgp)

* Share in JPY each year
egen isjypt = sum(jypt), by(year) // item sum
gen lnisjypt = ln(isjypt)

* generate country-wise value (JPY) shares
gen lnshp = lnjypt - lnisjypt // ln_share

* generate uv (unit value)
gen lnuvp = lnjypt - lnkgp // ln_price

gen d095=1 if year==1995
replace d095=0 if d095==.
gen d096=1 if year==1996
replace d096=0 if d096==.
gen d097=1 if year==1997
replace d097=0 if d097==.
gen d098=1 if year==1998
replace d098=0 if d098==.
gen d099=1 if year==1999
replace d099=0 if d099==.
gen d100=1 if year==2000
replace d100=0 if d100==.
gen d101=1 if year==2001
replace d101=0 if d101==.
gen d102=1 if year==2002
replace d102=0 if d102==.
gen d103=1 if year==2003
replace d103=0 if d103==.
gen d104=1 if year==2004
replace d104=0 if d104==.
gen d105=1 if year==2005
replace d105=0 if d105==.
gen d106=1 if year==2006
replace d106=0 if d106==.
gen d107=1 if year==2007
replace d107=0 if d107==.
gen d108=1 if year==2008
replace d108=0 if d108==.
gen d109=1 if year==2009
replace d109=0 if d109==.
gen d110=1 if year==2010
replace d110=0 if d110==.
gen d111=1 if year==2011
replace d111=0 if d111==.
gen d112=1 if year==2012
replace d112=0 if d112==.
gen d113=1 if year==2013
replace d113=0 if d113==.
gen d114=1 if year==2014
replace d114=0 if d114==.
gen d115=1 if year==2015
replace d115=0 if d115==.
gen d116=1 if year==2016
replace d116=0 if d116==.
gen d117=1 if year==2017
replace d117=0 if d117==.
gen d118=1 if year==2018
replace d118=0 if d118==.

gen lnxrall = ln(xrall)
gen lnxr = lnxrjpy - lnxrall

gen lnpppp = ln(pppp)
gen lncrop = ln(crop)
gen lnfamily = ln(family)
gen lnagex = ln(agex)

egen period = fill(1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 ///
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25)
gen lnperiod = ln(period)

xtset c_code year

*** ESTIMATION ***
xtivreg2 lnshp d095-d118 (lnuvp = l.key tr207), fe endog(lnuvp)
xtreg lnshp d095-d118 lnuvp, fe

nlcom 1-_b[lnuvp]
nlcom (F2018:exp(-((0)+_b[d118])/_b[lnuvp]))(F2017:exp(-((0)+_b[d117])/_b[lnuvp]))(F2016:exp(-((0)+_b[d116])/_b[lnuvp]))(F2015:exp(-((0)+_b[d115])/_b[lnuvp]))(F2014:exp(-((0)+_b[d114])/_b[lnuvp]))(F2013:exp(-((0)+_b[d113])/_b[lnuvp]))(F2012:exp(-((0)+_b[d112])/_b[lnuvp]))(F2011:exp(-((0)+_b[d111])/_b[lnuvp]))(F2010:exp(-((0)+_b[d110])/_b[lnuvp]))(F2009:exp(-((0)+_b[d109])/_b[lnuvp]))(F2008:exp(-((0)+_b[d108])/_b[lnuvp]))(F2007:exp(-((0)+_b[d107])/_b[lnuvp]))(F2006:exp(-((0)+_b[d106])/_b[lnuvp]))(F2005:exp(-((0)+_b[d105])/_b[lnuvp]))(F2004:exp(-((0)+_b[d104])/_b[lnuvp]))(F2003:exp(-((0)+_b[d103])/_b[lnuvp]))(F2002:exp(-((0)+_b[d102])/_b[lnuvp]))(F2001:exp(-((0)+_b[d101])/_b[lnuvp]))(F2000:exp(-((0)+_b[d100])/_b[lnuvp]))(F1999:exp(-((0)+_b[d099])/_b[lnuvp]))(F1998:exp(-((0)+_b[d098])/_b[lnuvp]))(F1997:exp(-((0)+_b[d097])/_b[lnuvp]))(F1996:exp(-((0)+_b[d096])/_b[lnuvp]))(F1995:exp(-((0)+_b[d095])/_b[lnuvp]))(F1994:exp(-((0))/_b[lnuvp]))

* Feenstra Implicit *****************************************
gen delta1 = lnuvp - lnshp*(.1596819)
xtivreg2 lnshp d095-d118 (lnuvp = delta1), fe endog(lnuvp)
*************************************************************
